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We study the spontaneous spin textures induced by magnetic dipole-dipole interaction in ferro¬ 
magnetic spinor condensates under various trap geometries. At the mean-field level, we show the 
system is dual to a Dirac string gas with a negative string tension in which the ground state spin 
texture can be easily determined. We find that three-dimensional condensates prefer a meron-like 
vortex texture, quasi one-dimensional condensates prefer the axially polarized flare texture, while 
condensates in quasi two dimensions exhibit either a meron texture or an in-plane polarized texture. 


I. INTRODUCTION 

The magnetic dipole-dipole interaction (DDI) is known 
to play a significant role in determining the long range 
behaviors of ultracold gases of spinful atoms, especially 
those with a large spin [lj-l3|. In the ferromagnetic spinor 
condensates of spinful bosons with nonzero spin expec¬ 
tation values (F) , the DDI usually leads to rich spin tex¬ 
tures that strongly depend on the trap geometry and the 
initial state preparation Meanwhile, the existence 

of spontaneous spin textures in dipolar spinor conden¬ 
sates has been verified by numerical simulations based 
on the mean-field Gross-Pitaevskii equation [§4l2j. How¬ 
ever, due to the complicated form of DDI, there still lacks 
an explicit and efficient theoretical way to understand the 
spin textures. 

In this letter, we develop an semi-analytical approach 
to the problem by introducing a duality at the mean-field 
level from a strong dipolar spinor condensate to a weakly 
interacting Dirac string gas with a negative string ten¬ 
sion. We find that the ground state is reached by forming 
as many closed strings with small enough curvatures as 
possible. Based on this principle, we are able to deter¬ 
mine the spin textures in various trap geometries used in 
cold atom experiments. We show the ground state spin 
texture of a large three-dimensional (3D) spherical con¬ 
densate is a meron-like vortex texture. Further, we show 
the axially polarized flare texture 0 is favored in quasi 
one-dimensional (ID) condensates, while in quasi two- 
dimensional (2D) pancake traps there is a phase transi¬ 
tion between an in-plane polarized texture and a meron 
texture driven by the pancake radius. In particular, mul¬ 
tiple merons can be seen to arise in quasi 2D traps when 
the aspect ratio is comparatively large, in agreement with 
the numerical results in Ref. [9j. The Dirac string gas 
picture thus offers a greatly simplified way to understand 
the spin textures in dipolar condensates. 


II. THEORETICAL FORMULATION 

We begin with the description of the ferromagnetic 
spinor condensates at low energies. In general, a spinor 
condensate of integral spin F atoms is characterized by 


a spinor order parameter of 2F + 1 components \H m (r) 
(to = —F, —F + 1, • • • ,F), which represents the coher¬ 
ent amplitude of annihilating a boson in the spin state 
| F,m) at position r jl3l - fl6l |. At the mean field level, 
the number of atoms per unit volume in the conden¬ 
sate is n(r) = J2m v ^™( r )^'m(r) = ’J , t(r)\E , (r), while the 
expectation of the local spin is (F(r)) = ip' (r)Ft^(r), 
where ip(r) = A>(r)/sjn(v) is the normalized spinor or¬ 
der parameter, and F is the spin matrix of the spin 
F representation. In a square well trap potential or 
in the central region of a harmonic trap within a size 
R <C n(r)/|Vn(r)|, we can approximate the particle 
number density n( r) = no &s & constant. Up to a 
U(l)xSU(2) transformation, the normalized spinor order 
parameter y>(r) is determined by the local interactions, 
which are characterized by F + 1 scattering lengths a 2 j 
between two atoms with total spin 2J (0 < J < F) 
0- A spinor condensate is said to be ferromagnetic if 
<p( r) = (0, • • • , 0,1, 0, • • • , 0, 0) in a certain basis, where 
the only nonzero component is to = Fq > 0. 

At low energies, the long wave length fluctuation of the 
order parameter gives rise to multiple gapless Goldstone 
modes. In a ferromagnetic spinor condensate, there is 
exactly one quadratic dispersion mode corresponding to 
the fluctuation of the local spin 0IH. In the presence 
of the magnetic dipole-dipole interaction, the effective 
low energy spin Hamiltonian of a ferromagnetic spinor 
condensate can be written as H = H 0 + Hjj, where 


H ° = JlI d3r ( VJP ( r )) 2 ’ 
H ° = \ 


3(^1 • ri2)(^2 • fi 2 ) 

4?r r\ 2 

(1) 


are the kinetic energy of the quadratic Goldstone mode 
and the DDI energy, respectively (see Appendix 0. We 
have used here the normalized local spin field F(v) = 
(F(r))/Fo that satisfies |.7 r (r)| = 1 for later conve¬ 
nience. M denotes the particle mass, a is defined by 
a = noh 2 [F(F + 1) — Fq]/ 4, and the interaction param¬ 
eter A is given by A = po(gFPB) 2 noFo, where gF is the 
Lande factor and ps is the Bohr magneton. In writing 
the Hd term, we have used the abbreviations fF, = Firf) 
and ri 2 = ri 2 ?i 2 = iq — r 2 , where fi 2 is the corresponding 
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FIG. 1. (a) A single column chain of head-to-tail spins is 

equivalent to a Dirac string with two monopoles ±Q m at the 
ends, (b) Under the dipolar interaction of spins, two open 
strings only have an interaction between the monopoles, while 
closed strings are free. 


unit vector. 

It is difficult to infer the low energy spin texture config¬ 
urations directly from the perspective of local spins, since 
the spin-spin interaction in Hp> is highly anisotropic. 
However, it is possible to understand Hp, more naturally 
in terms of Dirac strings. Consider a chain of atoms that 
are arranged uniformly along a curve as is shown in Fig. 
□I a ), with their spins aligned head-to-tail and parallel to 
the tangent of the curve. Since each atom carries a mag¬ 
netic dipole moment mi? = F), the curve can be 

exactly viewed as a Dirac string, with a positive magnetic 
monopole at one end and a negative one at the other. If 
we parameterize the curve as r(l) with an affine param¬ 
eter l , and require \dr/dl\ = 1 , the local spin field on 
the curve is then JF = dr/dl. We can easily calculate 
the dipolar interaction between such a Dirac string and 
a spin magnetic moment nip not belonging to the string: 


E mp — /To Qm J dl 

/ r f 


— / J’OQn 


4wr 3 47 rrf 


r m ff -3( m F-r)(r-f) 
47rr 5 


■ mi? = — mi? ■ B s 


( 2 ) 


where we have defined Q m = gp^sF^/h with the spacing 
between atoms on the string b ~ n 0 1 ^ 3 , while rj and 
r f are the positions of the ends of the string relative to 
m^. This means that effectively the magnetic moment 
mi? only feels the magnetic field of the monopoles at the 
ends of the string, which carry monopole charges ±Q m 
respectively. Remarkably, a closed Dirac string has no 
interaction energy with the other spins at all. 

Further, it is straightforward to prove that the interac¬ 
tion energy between two Dirac strings is solely given by 
the interaction between the four monopoles at the ends 
of the strings. Consider two strings L\ and L 2 as shown 


in Fig. QJb), whose ends are located at r 4 , 17 and r', r/ 
respectively. We shall parameterize string Li as r(/), and 
use the abbreviations r i(l) = r■ — r(/), r f(l) = r/ — r (l) 
to denote the displacement vectors from r (l) to the ends 
of string L 2 . The magnetic moment of an infinitesimal 
segment dl on string L\ is then mj? = gpg,BFodr(l)/b = 
Q m dv(l) = —Qmdri(l) = —Q m drf(l). According to the 
result of Eq. ©, the total DDI energy between the two 
strings is simply 


Eq(Li,L 2 ) — goQm 



r /( 0 

r *(0 

/ 

47rry(7) 3 

4nri(l) 3 


Mo Qm 


r ss r f (l)dr f (l) p p(Z)dr t (Q 
r if 47 T r f {l) 3 J rH 4nri(l) 3 


Mo Ql (J _1_ L + ±\ 

4 tt \r ff r if r fi r u ) 


( 3 ) 


where rff, rif , rfi and ru are the distances between 
monopoles at the ends of the two strings as shown in Fig. 
Hlb). The interaction energy between two monopoles is 
thus simply given by the Coulomb law. Accordingly, two 
parallel(anti-parallel) Dirac strings will repulse(attract) 
each other. Still, a closed string does not interact with 
any other strings. Therefore, the Dirac strings are almost 
free under DDI except for the monopole interactions be¬ 
tween their ends. 

In the continuum limit, we can view the flow lines 
of the local spin field lF(r) of the condensate as Dirac 
strings, each of which consists of a single column of 
atoms. We denote the cross sectional area of a string as 
da , which is around n 0 2 ^ 3 . The spacing between atoms 
on a string is then b = ( n^da ) _1 , and the monopole 
charge can be expressed as Q m = gF^BnoEoda, which 
is proportional to the cross sectional area da. 

So far we have seen the bulk of a Dirac string has no 
interaction with other strings at all under DDI. How¬ 
ever, as we shall shown in the below, a Dirac string has 
a negative self-interaction bulk energy coming from the 
DDI. To understand the origin of the self-interaction en¬ 
ergy, consider a straight Dirac string of length L. Since 
all the spins are aligned head-to-tail, there is an attrac¬ 
tive DDI energy — 2 (A/?rg)(l/ 47 rr 3 ) between any two of 
them, where r is the distance between them. This gives 
rise to the self-interaction energy of the string. In the 
limit b « L, the attraction energy felt by a spin on the 
string is approximately 



where we have used the fact 1/b 3 oc no- Thus, the self¬ 
attraction energy Ep 0 of each spin is a constant not sen¬ 
sitive to the string length L. 

The above estimation, however, cannot give the coeffi¬ 
cient in front of — A/no hi the expression of Ep 0 , since we 
do not know the exact relation between 1/b 3 and no in 
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the continuum limit. Instead, we can find out the coeffi¬ 
cient by directly computing the DDI energy of a simple 
concrete example of spin texture. Consider a 3D spheri¬ 
cal condensate with all spins polarized along z direction 
as shown in Fig. ©a). The total DDI energy E d felt 
by the spin at the center can be directly calculated in 
spherical coordinates from its definition: 


Appendix [Bj: 


Ei{Li,L 2 ) 


d(Tld(T2 


a 

M 




x 




r 2 ) , 


( 8 ) 


E d = — [ 2tt r 2 dr [ d6sm6- - ico f 6 = o , ( 4 ) 

n 0 Jo Jo 4ttH 

where R is the trap radius, and 9 is the polar angle. On 
the other hand, in the picture of Dirac strings, the DDI 
energy E d felt by the spin consists of the attraction en¬ 
ergy Ep 0 from its own string and the interaction with 
monopoles at the ends of all the strings derived in Eq. 
©. In cylindrical coordinates ( p,(j),z ), the cross sec¬ 
tional area of a string is da = pdpd(f> , and this energy is 
calculated as 

p\ rR / t>2 _ n 2 \ 

E d = E Fo + — 2npdp V P =E Fo + —. (5) 

n 0 J 0 47 tR 6 3n 0 

By equalizing Eq. (©) and Eq. ©, we see the attraction 
energy 


where is defined as V"j_ = V 2 — d and 14 = ri(?i), 
r 2 = Though looks complicated, this contact in¬ 

teraction simply implies that two nearby Dirac strings 
prefer to be coplanar, namely, \dri/dli — dr 2 /dl 2 |/|ri — 
r 21 —^ 0 as |r 1 — r 2 | —^ 0 . 

Therefore, we have seen the ferromagnetic dipolar con¬ 
densate at the mean-field level is dual to a Dirac string 
gas with a string self energy E str , a monopole interaction 
energy Eq and a contact interaction energy £ 7 , given in 
Eq. O, Eq. ([3]) and Eq. |5]), respectively. With E str 
mainly coming from the DDI Hp> and £/ coming from 
Hq, such a duality between individual spins and Dirac 
strings interchanges the roles of the kinetic energy and 
the interaction energy in some sense, and is thus a strong- 
weak duality. 


III. APPLICATION TO SPHERICAL TRAP 
GEOMETRIES 


The total self-interaction energy of a Dirac string of 
length L and of cross sectional area da is then 

E itl( L ) = -7T x n 0 Lda = ~XLda . ( 6 ) 

2 o 

This indicates that a Dirac string has a negative string 
tension T str = —Ada/6, and therefore the whole space 
will prefer to be filled with Dirac strings. The ac¬ 
tion for a single Dirac string can be written as S s t r = 
(' T str /h ) f dtdl = TstrE/h , where E is the worldsheet area 
of the string. The Dirac string therefore corresponds to 
a classical bosonic string theory with a U(l) gauge group 
01 - 

Now we turn to the kinetic energy Hq in Eq. m, and 
interpret it in the language of Dirac strings. By writ¬ 
ing the gradient operator into V = V|| + Vj_, where || 
and T stands for parallel and perpendicular to the string 
respectively, we can divide Hq into the following two 
parts: the first parallel part contributes an additional 
term to the string self energy. This is seen by noticing 
that V||^F = dT/dl = d 2 r/dl 2 , so the string self energy 
now reads 


f L „ 

1 , 

a 

( d 2 r A 2 ] 

/ dl 

—A + 

— 


1 0 

6 

M 

K dl 2 ) \ 


Note that \d 2 r/dl 2 \ is the curvature of the string. The 
second perpendicular part induces a nonnegative contact 
interaction between two Dirac strings L\ and L 2 (see 


The ground state spin texture can be more easily de¬ 
rived in the Dirac string gas picture. First, since an open 
string can always lower its monopole interaction energy 
by attaching its two ends together and forming a closed 
string, a strong dipolar condensate will prefer as many 
closed strings as possible. Further, to reduce the kinetic 
energy Hq, the close strings will tend to have a smaller 
average curvature and be locally coplanar to each other. 
These facts constitute the guiding rules for identifying 
the spin textures in a condensate under various trap ge¬ 
ometries. 

We first consider a dipolar condensate in a 3D spheri¬ 
cal square well trap potential of radius R. Based on the 
guiding rules above, the most natural spin texture one 
can expect is a meron-like vortex texture as shown in 
Fig. ©b). Since circular closed strings possess the small¬ 
est average curvature, the region p > Rq (in cylindrical 
coordinates) prefers to be filled with circular strings par¬ 
allel to each other, where Rq is a characteristic length 
scale to be determined. All the circular strings are sym¬ 
metric about the 2 axis, so that locally they are always 
coplanar. At p < Ro , the circular strings will have too 
large curvatures, and thus will gradually give way to spi¬ 
ral strings proceeding along the z direction to lower the 
kinetic energy Ho- Finally, the string becomes vertical 
and straight at the central axis p = 0. Consequently, 
there will be monopoles distributed on the p < Rq re¬ 
gions of the sphere’s surface, with a monopole area den¬ 
sity ua ~ EQm/da = EgppBnoFo. The total energy of 
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FIG. 2. (a) In the absence of dipolar interaction, all the spins 
will be polarized along the same direction, (b) The meron-like 
vortex spin texture of spherical dipolar condensates, (c), (d) 
and (e) shows several perturbation patterns that may arise in 
the spin texture, all of which increase the total energy. 
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FIG. 3. The ground state angle function £( p ) calculated for 
different values of Ro = (Sa/SXM) 1 ^ 3 R 1 ^ 3 . The dashed lines 
denote the positions of Ro- When Ro < R, Ro is a rather pre¬ 
cise boundary separating the C(p) ^ 0 region and the ((p) « 0 
region. 


the dipolar condensate can be estimated as 


spin texture ansatz, expressed in cylindrical coordinates 
(P, ft, z) as: 


E 3 D ~ 

6 


' p>Ro 


,0 a 1 

d 3 r—— + 

M p 2 


2/10(7^71,4) 2 
47ri?o 


A r 47ra , R 7r. 

-6 V+ ^r m ° g To + 2 XR > 


( 9 ) 


0 > 


where V = 47ri? 3 /3 is the volume of the condensate, 
and Rq/R is assumed small. The second term is the ki¬ 
netic energy Hq mainly coming from the curvatures of the 
strings, while the third term is the monopole interaction 
energy estimated by taking Rq as the average distance 
between monopoles [18]. By minimizing the total energy 
with respect to Rq, we get the characteristic length 

R 0 = (8a/3AM) 1 / 3 !? 1 / 3 . 


The condition for the meron-like spin texture to arise is 
then roughly R. > R 0l or R. > (8a/3AM) 1 / 2 . Therefore, 
the spin textures are more likely to arise for massive and 
large spin atoms. 

It is easy to see the meron-like spin texture is energet¬ 
ically stable. When perturbed, the circular strings may 
be tilted, misaligned, or distorted. In any case, the to¬ 
tal energy is raised higher. In the examples shown in 
Fig. He) and Fig. Hd), parallel open strings have to 
arise to keep the strings dense, which then cost an extra 
monopole interaction energy. In another case shown in 
Fig. [2] (e), no monopole occurs, but the strings are no 
longer locally coplanar to each other, which also increases 
the kinetic energy Hq. These perturbation patterns may 
have interesting dynamics, which is open to studies in 
the future. 

To verify the precision of our estimation above, we 
perform a numerical calculation based on the following 


•^"(r) = e <p cos </(p) + e z sin ((p) , (10) 

where £(/?) is a function of p only. This ansatz simply 
describes the spin texture shown in Fig. Hb). C (p) bas 
the physical meaning of the proceeding angle of a spiral 
string at radius p 1 and the string becomes circular when 
((p ) = 0. Fig. [3] shows C(p) calculated numerically for 
various values of Rq, and the details of the calculations 
can be found in Appendix[C] As is seen, our estimation of 
Ro gives a rather precise boundary between the region of 
spiral strings (£ (/>) ^ 0) and the region of circular strings 
(C(p) ~ 0)j even for Rq/R that is no longer very small. 


IV. APPLICATIONS TO CIGAR AND 
PANCAKE TRAP GEOMETRIES 

This analysis can be readily applied to the dipolar con¬ 
densates in other trap geometries. A practical and impor¬ 
tant geometry is the axially symmetric square well trap 
where the boundary is given by x 2 + y 2 + ( z/A) 2 = R 2 . 
By setting A 1, we get a quasi-ID cigar trap of 
length L c = 2 AR. In the extremely quasi-ID case, closed 
strings are no longer energetically favorable, since their 
curvatures are always too large. Instead, vertical strings 
along z direction are preferred, so that the amount of 
monopoles on the surface is minimal and the monopole 
energy is the smallest. In a realistic trap, the parti¬ 
cle number density no is not constant and usually lower 
near the surface. To validate our theoretical formulation 
in Sec. El we need to keep b = ( nodcr ) 1 constant for 
a string, so the effective cross sectional area density of 
strings n s ~ l/dcr oc no is also lower near the surface. In 
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FIG. 4. (a) The flare texture in a quasi-ID cigar trap, (b) 
The in-plane polarized texture in a quasi-2D pancake trap, 
(c) The meron texture in the pancake trap. 


the cigar trap case, the vertical strings are then forced 
to be bent outward to have a lower density near the sur¬ 
face, forming a flare texture as shown in Fig. Ha). By 
an energy estimation similar to what we did for spherical 
traps, one can show that circular closed strings begin to 
occur when A ~ (A M/a ) 1 ^ 3 or smaller, leading to 
a crossover from the flare texture to the 3D meron-like 
vortex texture. 

In the opposite limit A<Cl, the trap is a quasi-2D 
pancake of radius R. Correspondingly, there are two 
candidate ground state spin textures. The first spin tex¬ 
ture is shown in Fig. |4jb) , where all strings are in-plane 
and polarized along the same direction to minimize the 
amount of surface monopoles. The total energy of this 
texture consists of the string self-energy and the surface 
monopole interaction energy, or explicitly 

E& *~V + w\R 3 A 2 log j , (11) 

where V is the volume of the condensate, and w > 0 is 
a numerical factor. The second spin texture is the 2D 
meron texture shown in Fig. He). Similar to the 3D 
case, there is a characteristic radius R 2 D ■ The region 
outside R- 2 D is filled with in-plane circular strings, while 
inside R- 2 d strings become spiral and finally vertical at 
the center. In the case the half height of the pancake 
Zh = AR -C R 2 D 1 the total energy of the configuration 
can be estimated as 0 


pin) _ A t/ 

E 2D 


47TC* R 

-M^RTn 


— z h^2D 


( 12 ) 


Minimizing the energy we get R 2 D = (da/AM) 1 / 2 , which 
is independent of Zh or R. At a fixed half height Zh , one 
can show that the in-plane polarized texture is favored 
when R < i? c , where R c is a critical radius determined 
by solving E ^ = E while the meron texture arises 
when R > R c . We note the two textures are separated 
by a phase transition at R = R c instead of a crossover. 

When the pancake trap is not axially symmetric, mul¬ 
tiple merons may occur. As an example, in a rectangular 
trap with aspect ratio q = 2, the single-meron texture 



FIG. 5. In a large enough quasi-2D trap with an aspect ratio 
q = 2, the single meron in (a) costs more kinetic energy than 
the two antiparallel merons in (b). 


shown in Fig. H a ) costs a kinetic energy oc ZhR/R 2 D 
near the center. Instead, the two-meron texture in Fig. 
Hb) only costs a kinetic energy oc Zh \og(R/R 2 D), and 
thus will be favored at large R. The central spins of the 
two merons are aligned antiparallel to each other to fur¬ 
ther lower the DDI energy. Similarly, one may expect a 
g-meron spin texture to arise for a trap with aspect ratio 
q , which agrees well with the numerical results obtained 
in Ref. 


V. CONCLUSION 


In conclusion, we have shown that the duality to the 
Dirac string gas picture is an efficient semi-analytical way 
of determining the spontaneous spin texture in a dipolar 
spinor condensate. The total dipolar energy can be con¬ 
veniently estimated in the perspective of Dirac strings. 
Closed Dirac strings with curvatures as small as possible 
are generically preferred in the ground state. We expect 
this method to be further employed in future work to 
study the spin textures in external fields and the dynam¬ 
ics of spin textures. In addition, it will also be useful 
and intriguing to examine the possibility of constructing 
a quantum version of this duality as a generalization of 
the mean-field level duality presented here. 
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Appendix A: Mean field description of 
ferromagnetic spinor condensates 


Appendix B: The kinetic energy Ho rewritten in the 
Dirac string picture 


Spin F bosons with DDI are most generally described 
by a Hamiltonian H = H 0 + Hd , where Ho represents 
the kinetic energy and the local interaction [16j, and Hd 
is the long range DDI. In terms of the boson annihilation 
(creation) operators 0 = (0f, ■ ■ • , i/j-f) T , Ho is given by 


V 7 J =0 

(A!) 

where P 2 j = ® ip^)\2J,m){2J,m\(0 ® VO is the 

local two-particle interaction projected into the total 
spin 2 J channel, and g 2J is the corresponding interac¬ 
tion strength, and g 2J = 2nh 2 a2j /M in terms of the 
s-wave scattering length a 2 j. The whole Hamiltonian is 
17(1) x SU( 2) symmetric. In the mean-field approxima¬ 
tion, (if>m(r)) = ^m(r), where T m (r) = y/n(r)(p m (r) is 
the spinor order parameter. The mean field energy of Ho 
can then be written in the form 



Ho = err 


2 M 


|V \/n(r)ip(r)\ 2 + V{n(r),ip(r)) 


(A2) 

where V(n,ip) is invariant under U{ 1) x 51/(2) transfor¬ 
mations of ip. For a ferromagnetic spinor condensate, 
both the U(l) and SU(2) symmetries are spontaneously 
broken, and there are two Goldstone modes correspond¬ 
ing to the U(l) phase fluctuation and the spin direction 
fluctuation respectively [Til ]. At low energies, by inte¬ 
grating out the quantum fluctuation in density n(r) in 
the path integral formalism, we can reduce the mean- 
field energy (IA2I) to 



,3 r [ H 2 (d t (/>sF) 2 
[4 (d 2 V/dn 2 ) 
h 2 {dt<t>SF) 2 h 2 no 
4(d 2 V/dn 2 ) + 2M 


h 2 no 
2 M 


|V<p| s 


(V<V SF ) 2 + -(V^) 2 


(A3) 


where 0sf (r) is the U(l) phase, JF(r) = p^Fip/Fo is the 
normalized local spin, and the coefficient can be deter¬ 
mined to be a = noh 2 [F(F + 1) — Fo]/4. We note that 
the local interactions g 2J only contribute to the param¬ 
eter d 2 V/dn 2 in the first time dependent term, and do 
not enter the expression of the coefficient a. 

Similarly, the mean field DDI energy Hjj in Eq. © 
can be deduced by substituting the operator tp m with the 
mean field value 4/ m = y/nop m . 

Since the DDI only depends on the spin field r), 
the ground state always tends to have a constant U(l) 
phase 0s f provided the spin field ^-"(r) is smooth enough. 
Therefore, the effective low energy Hamiltonian of the 
system reduces to Eq. ©. 


To interpret the kinetic energy Hq in the Dirac string 
picture, we first note that it can be rewritten as 


Ho = ~/ d 3 r jr (r) . v 2 JP(r) 

= ~Jj J d 3 rid 3 r 2 S(r\ — r 2 ).F(r 2 ) • V 2 ^F(ri) , 


(Bl) 


where Vf is the Laplace operator with respect to iq. In 
the Dirac string picture, we can write the volume element 
d 3 rq as dlidcFi , and we have iF(iq) = diq/<iO, for i = 1,2. 
The kinetic energy then becomes 

H 0 = --J7 [ daida 2 dhdl 2 5 (v 1 -y 2 ) ( ^ . (B2) 
A 1 J dl 2 dti 

By partial integrating the expression two times, we find 
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dr 2 dr 1 


H 0 = - — / da 1 da 2 dlidl 2 -H L • —Vf<5(ri - r 2 
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dl2 dl 1 
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/ d(Jid(T2dlidl2 

1 ~dh 
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x [v^<j(ri - r 2 ) + df 1 S(r 1 - r 2 )] 
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= — / daida2dlidU 
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dr 2 diq 


x Vi<5(n -r 2 ) - — 
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dl 2 

daidli 
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= — / da\da2dl\dl2 
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1 - 


x V^_i5(ri - r 2 ) + 
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dh 
dr 1 

dh 

dr 2 dri 

~dh 

dl 2 r 


dl 2 

dadl 


dl 2 


d 3 ri 

Uf 


(B3) 


where Vj_ = V 2 — df. As is seen, the first term is a non¬ 
negative contact interaction between two Dirac strings, 
while the second term is an additional contribution to 
the self-energy of a Dirac string. 


Appendix C: Numerical calculations of spin texture 
in spherical traps 

In Eq. © we have derived an estimation of the to¬ 
tal energy of the spin texture in a spherical trap from 
the Dirac string perspective, and show that the meron 
texture has a core size Rq = (8a/3AM) 1 / 3 /? 1 / 3 . Here we 
present a more accurate numerical calculation of the con¬ 
figuration of the meron texture. Based on the reasoning 
given in the main text, the following spin texture ansatz 
is appropriate, written in spherical coordinates (p, 0, z) 
as: 

:F(r) = e^,cosC(p) +e z sin((p) , (Cl) 

bounded in the square well trap p 2 + z 2 < R 2 . As is 
shown in Fig. ©b), the monopole charges of this ansatz 
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are distributed on the surface of the trap, with an amount 
±n a pdpd(j> sin ((p) in the projected area pdpd(f>, where 
we have defined tia = gFpBnoEo■ The total energy of 
the spin texture is then the summation of the string self 
energy, the kinetic energy and the monopole energy: 


E^d = — TV - 
6 


a 


dT—(VFY + E Q 


=^fw^[( |r+ 


cos 1 2 3 4 5 6 7 8 9 10 £ 


+ / dpdp' 
Jo 


f 27T , Xpp' sin ( sin 

( * — TFT~ 


1 

~Dl 


1 

IT 


(C2) 


where 


D± = \Jr 2 - pp' cos <t> T Y(E 2 - P 2 )(-R 2 - P' 2 ) , (C3) 


and we have used the abbreviations £ = £(p), C = C(p')- 
\/2D± has the physical meaning of the distance between 
two monopole charges (+ for charges on the same hemi¬ 
sphere, — for charges on the opposite hemisphere). The 
integration j d(f>/D± can be expressed in terms of the 
elliptic K function to further simplify the energy func¬ 
tional. 

This energy functional of £(p) can be minimized via 
the steepest descent method. In practice, we start with 
an arbitrary function of £ (p), and add to the function an 
increment 


^C(p) — —eV(( P )E3D , 

repeatedly until the process converges, where e is a small 
number. Provided e is small enough, the numerical cal¬ 
culation converges well. The resulting functions £(p) for 
various values of Rq are shown in Fig. [3] 
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